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ABSTRACT 

We examine the possibility that very massive stars greatly exceeding the com- 
monly adopted stellar mass limit of 150Mq may be present in young star clusters in 
the local universe. We identify ten candidate clusters, some of which may host stars 
with masses up to 6OOM0 formed via runaway collisions. We estimate the probabilities 
of these very massive stars being in eclipsing binaries to be > 30%. Although most of 
these systems cannot be resolved at present, their transits can be detected at distances 
of 3 Mpc even under the contamination of the background cluster light, due to the 
large associated luminosities ~ lO^L© and mean transit depths of ^ lO^L©. Discov- 
ery of very massive eclipsing binaries would flag possible progenitors of pair-instability 
supernovae and intermediate-mass black holes. 
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1 INTRODUCTION 

Many observations support the statistical argument that 
the upper limi t to initial stellar masses is ~ 150M (7) for 
Pop II/I stars (|Figeill2005l : IZinnecker fc Yorkell2007l '). How- 
ever, this common n otion is challenged by t he recent spec- 
troscopic analyses of ICrowther et all (|201Cll ). in which star 
clusters NGC 3603 and R136 are found to host several 
stars with initial masses above this limit, including one star 
R136al with a current mass of ~ 265Mq. Also, candidate 
pair-instability supernovae, which require progenitors with 
masses above 200M(7), have been o bserved in the low red- 
shift universe l|Gal-Yam et al.ll2009l ). Therefore, it is worth 
exploring methods to confirm the existence of a very massive 
star (VMS), defined here as a star with a stellar mass signif- 
icantly greater than the stellar mass limit, i.e. M > 2OOA/0. 

Unless the VMS is very close by, it is extremely 
difficult to spatially resolve the VMS from stars in its 
vicinity. Indeed, the central component of R136 was once 
though t to be an extremely m assive ^, IO^Mq star 
Cassinelli. Mathis fc avageHToill ). efore IWeigelt fc Baieil 



19851 ) resolved it as a dense star cluster via speckle inter- 



ferometry. As for spectroscopic measurements, verification 
of a single VMS is further complicated by the fact that the 
effective temperature Te// of Pop I stars above 10^ Mq de- 
pends very weakly on mass, wit h log(re///K) « 4.7-4.8 
iBromm. Kudritzki fc Loebl I2OO1I ) for stars between 10 - 
10^ Mq. Moreover, a hot evolved star with an initial mass 
below 10^ Mq can nevertheless reach these temperatures in 
its post main-sequence evolution and mimic a VMS. 



The most accurate method of constraining the stellar 
masses of distant stars is by measuring the radia l velocity 
and l i ght curves of the star in an eclipsin g binary (|Bonanosl 
I2OO9I : iTorres. Andersen fc Gimenej I2OIOI '). The light curve 
provides a wealth of information about the binary, includ- 
ing its orbital period, inclination, eccentricity, as well as the 
fractional radii and flux ratio of the binary members. The 
radial velocities found from a double-lined spectroscopic bi- 
nary further provide the mass ratio of the binary. With the 
above information, the individual masses of each star in the 
binary can be calculated via Kepler's third law. Searches for 
massive eclipsing binaries in s tar clusters within our ow n 
Galaxy are already underway (|Koumpia fc BonanosI [2OI1I ) . 
and techniques have been sug gested for binary searches in 
other galaxies (|Bonanosll2"0ia) . 

In this Letter, we estimate the masses and properties 
of VMSs that may have formed via collision runaways in 
a number of very young, dense, and massive star clusters 
in the local universe. We calculate the probability of these 
VMSs to be in eclipsing binaries, and find their expected 
transit depths and observability. 



2 VERY MASSIVE STARS 

Shortly after a dense star cluster forms, its most mas- 
sive constituents sink to the center via dynamic friction 
and form a central subsystem of massive stars. In suffi- 
ciently dense environments, these massive stars may un- 
dergo runaway collisions and merge into a single VMS 
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llGiirkan. Freitag fc Rasid l2004l : iFreitag. Giirkan fc Rasic 



to 



WMq. iPortegies Zwart et al 



l2006f) . possibly up 
(|200a ) gives a fitting formula for the stellar mass rrir of the 
final runaway product, calibrated by N-body simulations for 
Salpeter-like mass functions: 



0.01 Mc 1 + 



lOOMyr 



where t^h is the relaxation time, 



irh 



f Mc 



(m) 



(1) 



(2) 



Here Mc is the cluster mass, r^ir is its virial radius, and 
(m) « O.5M0 is the average stellar mass. 

Using the compilation of stars clusters in 
the local universe and the ir prop erties from 
IPortegies Zwart. McMillan fc Giele^ (|2010l ). we have 
listed in Table [T] several young, dense star clusters that 
may host a runaway collision product of mass > 2OOM0 
which may have not yet ended its life as a star. We restrict 
our sample to clusters with mean determined ages younger 
than 3.5 Myr. This may already be insufficiently selective, 
as stars born with masses > 200Mp, are expect ed to have 
lifetimes of only 2-3 Myr (|Yungelson et al]|2008l ): however, 
in the runaway collision scenario, the VMS builds up 
its extraordinary mass via mergers over ~ 1 — 2 Myr, 
and therefore its host cluster may have an age exceed- 
ing the 2-3 Myr limit. Of course, these observed cluster 
properties should not be taken as certain; for example, 
lUbeda. Mafz-Apellaniz fc MacKentvl (2007) find the ages of 
NGC 4214 I-A and I-B to be ~ 4-5 Myr, likely too old for 
a VMS to be present. Conversely, there may be candidate 
clusters with VMSs that we have missed. The predicted 
runaway masses are only approximate, but give a sense of 
the mass range of VMSs that may lurk at the center of 
these very young and dense clusters. 

Alternatively, if feedback effects are moderate, it may 
be possible for a protostar to grow without a fixed mass limit 
via mergers or via the accretion of extremely dense gas. In 
this case, the mass of the most massive star m„ formed in a 
molecular cloud scales with the mass of that cloud, and thus 
will be correlated with the mass of its eventual host cl ust er 
l|Larsonlll983 . l2003l : IWeidner. Kroupa fc Bonnelj|20ld ): 



1.2 A/c 



(3) 



If the above relationship is valid for cluster masses > 5 x 
1O*M0, VMSs will not be restricted to dense clusters, since a 
collision runaway is no longer necessary for achieving masses 
> I5OM0 (see Table [T]). 



3 ECLIPSE PROBABILITY 

The fraction of massive O-type sta rs in binaries /;, is ob- 
served to be extremely high > 70% JChini et al. 2012 ), and 
approaches 100% in some environm ents (jMason et al TI2OO9I : 
iBosch. Terlevich fc Terlevichll2009l '). Although there is no re- 
lated observational data on VMSs, numerical simulations 
indicate that the collision runaway product in young, dense 
star clusters is generally accompanied by a companion star 
(Portegies Zwart, private communication). 



As the period distribution for our hypothetical VMS bi- 
naries is unknown, we assume their periods share the same 
cumulative distribution function (CDF) as the periods of 
massive binaries determined from observations. The CDF of 
the orbital period (p, in days) for massive binaries follows 
a 'broken' Opik law, i.e. a bi-uniform distrib ution in logp, 
with the break at p = 10 (|Sana fc Evandlioi lV There is an 
overabundance of short period binaries, with 50% to 60% of 
binaries having periods less than 10 days. The correspond- 
ing probability distribution function PDF{p) of the orbital 
period is: 



A, for 10°-^ p s; 10 



PDF(p) 



In 10 



5p ' 



for 10 < p ^ lO'' 



with the normalization J PDF{p)dp = 1. 

By integrating over uniformly distributed inclinations, 
it is easy to show that the eclipsing probability of a binary 
system at any depth is Pe(a) — -St, where Rt — R\ + R2 is 
the sum of the radii of both components in the binary, and 
a is the orbital distance. From Kepler's third law, we can 
express the eclipsing probability as a function of p instead: 



Pe{p) = Rt 



f)',GM. 



(■1) 



where Mt = Mi + M2 is the total system mass. Therefore, 
integrating over the period distribution, the probability that 
a massive binary will be an eclipsing binary to an observer 
on Earth is 



Pe(p) PDF{p) dp 



0.053 



Rt 



Rg 



Mt 



Me- 



(5) 



For convenience, we ignore any effects of eccentricity; tidal 
evolution will rapidly circularize the orbit for binaries with 
periods below p = 10 days, which account for 88% of the 
above eclipsing systems. Dynamical effects would harden 
a wide-separation massive binary system in the core of a 
dense cluster on a timescale much shorter than 1 Myr. Since 
three-body interactions tend to eject the lightest star, the 
companion to the VMS will likely be a massive star, though 
not as massive as the runaway product. 

The large radii of VMSs coupled with their high binary 
fraction (and short period binaries being common), imply 
significant eclipsing probabilities for VMSs. Using R136al 
as an example of the primary star, with a radius ~ 35-R0 , 
and a secondary Sun-like star, the eclipsing probability is 
29%, while for a more massive secondary star more common 
in the core of a young massive star cluster, e.g. a BO star of 
mass ~ I8M0 and radius ~ 7Rq, the eclipse probability is 
34%. Note that the eclipsing binary probability in equation 
(O is not sensitive to the secondary star parameters, as long 
as its radius is small relative to the primary. 

Assuming a companion BO star, we list the eclipsing 
binary probabilities for our candidate VMSs in Table [TJ cal- 
culated from equation ((5)1, except that we limit the inte- 
gration over p to periods corresponding to orbital distances 
exceeding both the radius of the VMS and the Roche limit 
for the companion. This restriction reduces Pe, and leads 
to the larger VMSs having slightly smaller eclipsing proba- 
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Table 1. Possible very massive stars in star clusters and their eclipse probabilities. The predicted runaway collision product mass 
rrir is calculated from equation {T}. Another possible VMS stellar mass rrtu is calculated via the relationship between the cluster 
mass and its most massive star in equation JS}. All masses are in units of Mq, the cluster age is measured in Myr, and the virial 
radius ryir is in units of pc. If we optimistically choose the largest m ass of and m„ for the pr i mary mass Mi, we can calculate 
its luminosity Li (in Lq) and radius _Ri (in Rq) using the models of lBromm. Kudritzki fc Loebl 1I2OOII ). assuming a characteristic 
stellar metallicity (Z/Zq) = 0.3. We calculate the eclipsing probability Pg assuming that the companion is a BO star, although 
the result is weakly sensitive to the companion mass. For generality, the expected transit depth (S) is averaged over a uniform 
distribution in the binary mass ratio q, up to a companion mass of W'^Mq, assuming non-grazing orbits, i.e. <5 (i?2/Ri)^. For all 
VMS candidates below, the expected dip in luminosity from the eclipse is ~ IO^Lq. 



Galaxy Name Ref Age logMc r^ir rrir ruu Li Ri Pe (S) 



Milky Way 


Arches 


1 


2.0 


4.30 


0.68 


192 


103 


5e6 


44 


39% 


16% 


LMC 


R136 


2,3,4 


3.0 


4.78 


2.89 


406 


170 


le7 


61 


36% 


8% 


SMC 


NGC 346 


5 


3.0 


5.60 


15.28 


640 


397 


2e7 


76 


34% 


5% 


M33 


NGC 604 


6 


3.5 


5.00 


48.21 


97 


213 


6e6 


46 


38% 


15% 


NGC 1569 


C 


6 


3.0 


5.16 


4.50 


672 


252 


2c7 


77 


34% 


5% 


NGC 4214 


I-A 


6 


3.5 


5.44 


28.69 


305 


337 


le7 


56 


36% 


10% 


NGC 4214 


I-B 


6 


3.5 


5.40 


9.85 


619 


323 


2e7 


74 


34% 


6% 


NGC 4214 


II-C 


6 


2.0 


4.86 


23.43 


129 


185 


5e6 


43 


39% 


17% 


NGC 4449 


N-2 


6 


3.0 


5.00 


3.57 


565 


213 


2e7 


71 


35% 


6% 


NGC 5253 


IV 


6 


3.5 


4.72 


5.26 


271 


160 


8e6 


51 


37% 


12% 



fP lFiger. McLean fc Morrij l|l999l l: (2} 'Hunter et al.' ('1995"); (•3VMackev fc Gilmorj l|2003h : (41 1 Andersen et al.l l|200g^ : (5) 

ISabbi ct al. (2008); (61 .Mai'z-Apellaniz 12001^ . 



bilities; nevertheless, the echpsing probabilities for all VMS 
candidates exceed 1/3. 



4 OBSERVABILITY OF TRANSIT 

VMSs have spectacular luminosities in the range of Lq; 
for example, R136al is observed to have ~ 8.7 x 10*' Lq. 
Even at a distance of 3 Mpc - roughly the distance of the 
farthest host galaxy in Table [1]- a star like R136al would 
still have an apparent bolometric magnitude of 14.8. How- 
ever, VMSs with Teff ~ 5 X 10"* K emit primarily in the 
ultraviolet, requiring bolometric corrections of BC ~ 4.6. 
Still, such a VMS will be within the V-band limiting magni- 
tude of ground-based 1-meter telescopes. For the VMS can- 
didates in Table [1] with a hypothetical BO-star companion, 
the transit depth exceeds 10^ in all cases, which at 3 Mpc 
is just within the single-visit limiting magnitudes of future 
synoptic surveys such as Pan- STARR^I and the Large Syn- 
optic Survey Telescopfl Of course, given the shortlist of 
host clusters in Table [U one can use deep, targeted obser- 
vations of the individuals clusters with existing telescopes, 
instead of uniform field surveys. 

However, in massive binaries, the mass ratio between 
the primary and secondar y star q = Mjj/Mi is observed 
to have a flat distribution (|Sana fc EvansI 201 lI V Unlike the 
transit probability, the transit depth is very sensitive to the 
companion star radii, so using a BO star as the companion 
may be overly conservative. Since only one VMS is expected 
to form in the collision runaway scenario, here we assume the 
distribution of companion star masses is uniform between 
1 to 100 Mq. Using typical mass-radius relationships, we 
show in Table [T] the expected transit depth {5} integrated 
over the range of companion star radii. Figure [T] illustrates 



^ http: / /pan- Starrs . if a.haw aii ■ edu/public/| 
^ http : //www . Isst . or g/lsst / I 



sample light curves for a VMS binary at 3 Mpc with different 
companion star masses and radii at different inclinations. 

For clusters outside the Milky Way and its satellites, it 
is currently impossible to resolve a VMS from other massive 
stars in a dense cluster core. Hence, we consider the luminos- 
ity of the host cluster as a contaminating third hght source 
to the eclipsing binary hght curve. If the VMS is present, 
it will contribute a significant fraction of the bolometric lu- 
minosity of the cluster (at least 10% and exceeding 50% in 
some cases), and an even larger fraction of the UV flux. The 
integration time t needed to reach a target signal-to-noise 
ratio SNR for detecting a transit can be approximated as: 



t 



6 seconds x 

fband 



c 



0.2 



108 L, 

1 r 







d 



10 eV 



3 Mpc 
A 



4 X 104 cm2 



SNR 


2 


fvMS 


-2 


■ 5 ' 


10 




0.1 




10% 



(6) 



where Lc is the bolometric luminosity of the cluster, d is the 
distance to the cluster, fband is the fraction of total flux that 
is observed (due to the spectral energy distribution, filter 
bandpass, CCD response, atmospheric transmission etc.), 
Eband is the characteristic observed photon energy, A is the 
collecting area of the telescope, fvMS is the fraction of total 
observed flux from the VMS primary, and again S is the 
transit depth. 

Note that the Hubble Space Telescope (HST) would 
collect > 10" UV photons per second from a lO^L© VMS 
even at a distance of 3 Mpc, thus detecting a J ~ 10% transit 
depth at SNR = 10 in tens of seconds of integration time. 
Obscuration by dust along the line-of-sight may reduce the 
observed UV flux. For V-band observations, a very young 
~ 10^ Mq cluster can be as bright as My ~ —12, while a 
3OOM0 VMS will have Mv ^ -8, i.e. the VMS wiU only 
contribute fvMS ~ 2.5% of the cluster hght in the visible 
band. Nevertheless, a 2-meter ground-based telescope will 
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Figure 1. Example light curves for a VMS eclipsing binary. The 
primary has parameters similar to R136al, while the secondary is 
either a 18Mq (dashed line) or IOOM0 (straight line) star, with 
appropriate radii and luminosities. The apparent magnitude m 
(bolometric) is plotted for these systems at 3 Mpc. The thick and 
thin lines correspond to inclinations of 90° and 70° , respectively; 
the period is 5 days in both cases. Reflections and limb - darke ning 
using the model of iDiaz-Cordoves. Claret &: Gimene are 
taken into account, but ellipsoidal variation is ignored. 



need less than an hour of integration time to detect the 
transit, which is eminently feasible as the transit duration 
r ^ p{Ri/-Ka) (X p^^^ for a VMS eclipsing binary will be 
> 10 hours for all relevant orbital periods. 

Other less massive eclipsing binaries in the host clus- 
ter will also contaminate the hght curve, but their transit 
depths will likely be negligible compared to the VMS's lumi- 
nosity. Non-binary random occultations of the central VMS 
can replicate a large transit depth, but using a King model 
for the cluster density profile (King 1966), we find these 
events occur less than once every 10*' years. 



5 STELLAR MASS DETERMINATION 

The extraordinary luminosity of a VMS should allow its ra- 
dial velocity to be measured. However, the mass ratio g, 
critical for model-independent determination of the individ- 
ual masses, can only be found when the radial velocities are 
determined for both components of the binary. Such double- 
lined spectroscopic binaries are easily observable when the 
components h ave similar luminosities, within a factor of 5 
of each other (jKallrath fc Milonell2009l ). As the luminosity 
of massive stars near the Eddington limit scales with mass, 
this criteria roughly corresponds to q > 0.2, which for an 
uniform distribution in g £ (0, 1) is quite likely to occur. 

Nevertheless, if the companion is small, and only spec- 
tral lines from the VMS are detected, then the mass ratio 
q cannot be unambiguously obtained. Instead, the mass of 
the VMS can be expressed as a single function of q: 

Mi= ^^+3'^^' .\j {M,,M2,i), (7) 
q'^ sm t 

where f{Mi,M2,i) is the mass function, which can be cal- 
culated using quantities derivable from the spectroscopy of 



a single-lined spectroscopic binary, and the inclination i is 
derivable from the eclipsing binary Ught curve. Unfortu- 
nately, equation varies sharply as oc q~^ for g <g 1. Since 
g can be as small as ~ 0.01 for VMSs in Table [T] crude con- 
straints on the mass ratio, e.g. g < 0.2 (when light from the 
secondary is not observed) cannot establish tight minimum 
stellar mass constraints on the VMS primary. 

However, since a total eclipse S — >■ 100% is extremely 
unlikely given the large radii of VMSs, if the mean value 
~ IO^Lq dip in the hght curve is in fact observed, it will 
immediately imply the existence of a star > 10^ M©. Hence, 
although sophisticated light-curve fitting with stellar models 
would be required, eclipsing single-lined spectroscopic bina- 
ries still offer an attractive avenue for inferring the presence 
of a VMS greatly exceeding the 150 Mq stellar mass limit. 



6 DISCUSSION 

A search for periodic fiux variations (as shown in Fig. 1) 
due to transits of the VMS candid ates in Table [T| would 
be of considerable interest. Although ICrowther et all(|2010l ) 
made robust arguments against R136al being a wide sepa- 
ration binary or an equal-mass binary, this source could still 
involve a short-period, unequal-mass binary system. The 
Arches cluster is observed to have no star s curre ntly above 
the I5OM0 mass limit, but lCrowther et al.l (|2010l ) also found 
with contemporary stellar and photometric results that the 
most luminous stars in the Arches cluster had initial masses 
approaching 2OOM0. 

The radii of VMSs are depende nt on their metallici- 
ties and rotation (|Langer et al.lT2007l ). If the VMS radii in 
Table [T] were smaller by ~ 25% (e.g. at much lower metal- 
licities), all listed eclipsing probabilities would still remain 
above 1/3, but the expected transit depth would increase up 
to (S) ^ 20%. As for the companion star, for most O stars, 
the point of unity Thomson optical depth occurs close to 
the hydrostatic radius, but when stellar mass loss exceeds 
~ 1O~*M0 yr~^, the photosphere r ~ 1 occurs in the wind 
itself, effectively increasin g the star's radius. This o ccurs for 
Wolf-Rayet companion s (ILamontagne et al. l '1996') and for 
companions > 6OM0 (|Vink. de Koter fc L amcrs 2000), in 
which case our eclipse probabilities and transit depths are 
too conservative. 

Binaries can be broadly classified into detached systems, 
where neither component fills its Roche lobe, versus semi- 
detached or over-contact systems, where at least one com- 
ponent exceeds its Roche lobe. VMSs in detached binaries 
have much more sharply defined eclipses, and more impor- 
tantly, they do not undergo mass transfer and lose mass to 
their companion. To find the probability that our VMS can- 
didates in Table [T] are detached eclipsing binaries, we limit 
the integration in equation ((5]) to periods p > 5 days, cor- 
responding to orbital distances where t he Roche lobe o f the 
VMS is always greater than its radius (|Eggletonlll983 ). For 
our VMS candidates, the detached eclipsing binary proba- 
bility is ~ 17%, i.e. roughly half of all eclipsing systems. 

However, non-pristine massive stars can also lose mass 
via strong winds driven by radiation pressure, with a mass 
loss rate increasing with metallicity. Post main-sequence 
VMSs can also lose mass eruptively or via pulsational in- 
stabilities, although mass loss near the end of the star's life 



© 2012 RAS, MNRAS 000, [THT] 



Very massive stars in eclipsing binaries 5 



(e.g. the pulsational pair-instability) is not likely to change 
the observability of our VMS candid ates. Under extraordi- 
nary mass loss via winds, ,Glcbbcck et al.l (|2009l ) found the 
highest mass attained by a collision runaway product to be 
~ 400Mq , although the star remaine d at this mass range 
for only ~ 0.2 Myr. On the contrary, ISuzuki et al] ()2007l l 
found that stellar mass loss does not inhibit the formation 
of a VMS of ~ W^Mq. 

If VMSs do in fact form via collision runaways in young, 
dense star clusters, and retain sufficient masses at the end 
of their lives, t hey may explode as p air-instability super- 
novae (PISNe) l|Yungelson et al.l l2008Y The creation rate of 
runaway products is in fact consistent with the current ob- 
served PISN rate (|Pan. Loeb fc Kasenll2012l '). However, the 
most massive VMSs may collapse directly into an interme- 
diate mass black hole (IMBH) via the pho todisintegration 
instability (|Wooslev. Heger fc Weaveil 120021 '). Tentative ev- 
idence has been claimed for IMB Hs at the center of old 
globular clusters (|Lou fc Wull2012l '). and extragalactic ultra- 
luminous x-rays sources associated with young star clusters 
l|Ebisuzaki et al.l |2001; Farrc U ct al. 2009). The identifica- 
tion of VMSs that can serve as the progenitors of PISNe 
and IMBHs will help move these extreme astrophysical ob- 
jects from the realm of speculation into reality. 



ACKNOWLEDGMENTS. 

We thank Paul Crowther, Dave Latham, Philip Myers, 
Guillermo Torres, and Simon Portegies Zwart for helpful 
discussions. We thank Eran Ofek for his eclipsing binary 
script^, which we built upon to generate our example light 
curves. TP was supported by the Hertz Foundation. This 
work was supported in part by NSF grant AST-0907890 and 
NASA grants NNX08AL43G and NNA09DB30A. 



REFERENCES 

Andersen M., Zinnecker H., Moneti A., McCaughrean 
M. J., Brandl B., Brandner W., Meylan G., Hunter D., 
2009, ApJ, 707, 1347 

Bonanos A. Z., 2009, ApJ, 691, 407 

— , 2012, in lAU Symposium, Vol. 282, lAU Symposium, 
pp. 27-32 

Bosch G., Terlevich E., Terlevich R., 2009, AJ, 137, 3437 
Bromm V., Kudritzki R. P., Loeb A., 2001, ApJ, 552, 464 
Cassinelli J. P., Mathis J. S., Savage B. D., 1981, Science, 
212, 1497 

Chini R., Hoffmeister V. H., Nasseri A., Stahl O., Zinnecker 

H., 2012, ArXiv e-prints 
Crowther P. A., Schnurr O., Hirschi R., Yusof N., Parker 

R. J., Goodwin S. P., Kassim H. A., 2010, MNRAS, 408, 

731 

Diaz-Cordoves J., Claret A., Gimenez A., 1995, A&AS, 110, 
329 

Ebisuzaki T. et al., 2001, ApJ, 562, L19 
Eggleton P. P., 1983, ApJ, 268, 368 

Farrell S. A., Webb N. A., Barret D., Godet O., Rodrigues 
J. M., 2009, Nature, 460, 73 

^ |http : //«lse-obs ■ tau ■ ac ■ 11/- eran/matlab . html | 



Figer D. F., 2005, Nature, 434, 192 

Figer D. F., McLean I. S., Morris M., 1999, ApJ, 514, 202 
Freitag M., Giirkan M. A., Rasio F. A., 2006, MNRAS, 368, 
141 

Gal- Yam A. et al., 2009, Nature, 462, 624 

Glebbeek E., Gaburov E., de Mink S. E., Pols O. R., Porte- 
gies Zwart S. F., 2009, A&A, 497, 255 

Giirkan M. A., Freitag M., Rasio F. A., 2004, ApJ, 604, 
632 

Hunter D. A., Shaya E. J., Holtzman J. A., Light R. M., 
O'Neil, Jr. E. J., Lynds R., 1995, ApJ, 448, 179 

KaUrath J., Milone E. F., 2009, Eclipsing Binary Stars: 
Modeling and Analysis 

King I. R., 1966, AJ, 71, 64 

Koumpia E., Bonanos A. Z., 2011, in lAU Symposium, Vol. 

272, lAU Symposium, Neiner C, Wade G., Meynet G., 

Peters G., eds., pp. 515-516 
Lamontagne R., Moffat A. F. J., Drissen L., Robert C, 

Matthews J. M., 1996, AJ, 112, 2227 
Langer N., Norman C. A., de Koter A., Vink J. S., Cantiello 

M., Yoon S.-C, 2007, A&A, 475, L19 
Larson R. B., 1982, MNRAS, 200, 159 
— , 2003, Reports on Progress in Physics, 66, 1651 
Lou Y.-Q., Wu Y., 2012, MNRAS, 422, L28 
Mackey A. D., Gilmore G. F., 2003, MNRAS, 338, 85 
Mai'z-ApeUaniz J., 2001, ApJ, 563, 151 
Mason B. D., Hartkopf W. I., Gies D. R., Henry T. J., 

Helsel J. W., 2009, AJ, 137, 3358 
Pan T., Loeb A., Kasen D., 2012, MNRAS, 2979 
Portegies Zwart S. F., Baumgardt H., McMiUan S. L. W., 

Makino J., Hut P., Ebisuzaki T., 2006, ApJ, 641, 319 
Portegies Zwart S. F., McMillan S. L. W., Gieles M., 2010, 

ARA&A, 48, 431 
Sabbi E. et al., 2008, AJ, 135, 173 

Sana H., Evans C. J., 2011, in lAU Symposium, Vol. 272, 

lAU Symposium, C. Neiner, G. Wade, G. Meynet, & 

G. Peters, ed., pp. 474-485 
Suzuki T. K., Nakasato N., Baumgardt H., Ibukiyama A., 

Makino J., Ebisuzaki T., 2007, ApJ, 668, 435 
Torres G., Andersen J., Gimenez A., 2010, A&A Rev., 18, 

67 

Ubeda L., Mafz-ApeUaniz J., MacKenty J. W., 2007, AJ, 
133, 932 

Vink J. S., de Koter A., Lamers H. J. G. L. M., 2000, A&A, 
362, 295 

Weidner C, Kroupa P., BonneU I. A. D., 2010, MNRAS, 
401, 275 

Weigelt G., Baier G., 1985, A&A, 150, L18 

Woosley S. E., Heger A., Weaver T. A., 2002, Reviews of 
Modern Physics, 74, 1015 

Yungelson L. R., van den Heuvel E. P. J., Vink J. S., Porte- 
gies Zwart S. F., de Koter A., 2008, A&A, 477, 223 

Zinnecker H., Yorke H. W., 2007, ARA&A, 45, 481 



© 2012 RAS, MNRAS 000, [iHI] 



